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Theoretical predictions underlying determinations of the fine structure constant a and the electron-to-proton 
mass ratio m e /m p are reviewed, with the emphasis on the bound electron magnetic anomaly g — 2. The theory 
of the interaction of hydrogen- like ions with a magnetic field is discussed. The status of efforts aimed at the 
determination of O (a(Za) 5 ) and O (a 2 (Za) 5 ) corrections to the g factor is presented. The reevaluation of 
analogous corrections to the Lamb shift and the hyperfine splitting is summarized. 



1. Introduction 

The world of our everyday experience is largely 
controlled by two dimensionless fundamental 
quantities: the fine structure constant a ~ 1/137 
and the electron-to-proton mass ratio m e /m p ~ 
1/1836. The former describes the strength of the 
electromagnetic interaction. The two together re- 
late the energy of atomic excitations to the bulk 
of the atomic mass. 

Among the myriads of phenomena governed by 
these two constants, only a handful are suitable 
for their precise determination. Two character- 
istics are required: they have to be amenable 
to measurements and they must be theoretically 
very well understood. Most of such systems con- 
sist of only one, two, or three bodies. 

In this talk we focus on the theoretical aspects 
of hydrogen-like ions. Measurements of their in- 
teraction with the magnetic field are a crucial in- 
gredient in the best determination of m e /m p , and 
indirectly help to determine a with the second- 
best accuracy. 

We first briefly review the present knowledge of 
both constants, to illustrate the level of accuracy 
already achieved. A system of ultimate simplicity 
is a single electron. It is studied with the Penning 
trap where it forms what is sometimes called a 
geonium atom |T| . The spectrum of lowest lying 
states is sensitive to the electron anomalous mag- 
netic moment, g — 2, and thus to a through the 



theoretical expression 






Ac, 



(1) 



where C\ are the coefficients determined from an 
i-loop correction to the interaction of an electron 
with a weak homogeneous magnetic field, calcu- 
lated in quantum electrodynamics (QED); and 
Aa are small corrections due to weak interactions 
and hadronic effects. The QED effects are now 
being studied at five-loops |2l3j . Indeed, the miss- 
ing information about C5 is the leading error in 
the present best value of a [3] , 



a" 1 (g - 2) = 137.035 999 084(39)th(33) e 



(2) 



The relative uncertainty of this evaluation is 0.37 
part per billion (ppb). 

The Penning trap can also be used to measure 
m e /m p by comparing the cyclotron frequency 
lo c = eB/m of an electron with that of a proton. 
The best value obtained in this way is accurate 
to 2.2 ppb [5], 



m. 



(free e") = 1836.152 666 5(40). 



(3) 



The systematic limitation of this determination 
is the relativistic shift of the electron mass due to 
thermal excitations of the magnetron motion of 
the electron. For this reason, studies of electrons 
bound in a heavier ion, that will be discussed be- 
low, are superior: the thermal velocity of such a 
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heavy system is eliminated as the source of rele- 
vant errors. 

The second-best method for finding a is based 
on the excellent knowledge (to within 0.007 ppb) 
of the Rydberg constant [6] , 

o 

Roo = e = 10 973 731.568 527(73)/m, (4) 

where h denotes the Planck constant and c the ve- 
locity of light. The latter being an exact number, 
a can be determined if the ratio of the electron 
mass to the Planck constant is measured. Such a 
direct measurement is difficult, but it is possible 
to measure the ratio of an atom mass m a t om to 
h for atoms such as cesium and rubidium [7 8 . 
Then, if m e /m p and TOp/m a t om are known, one 
can determine a. (This underscores again the im- 
portance of m e /m p .) The best value obtained in 
this way, with rubidium atoms, is 



aT 1 (Rb) = 137.035 999 450(620), 



(5) 



a 4.6 ppb determination. Efforts to reduce this 
error to about 1 ppb are being undertaken [9]. 

A macroscopic effect that competes with the 
atomic scale phenomena in the possible accuracy 
of determining a is the quantum Hall effect |10j , 
giving the value 



a -1 (QHE) = 137.036 003 700(3300), 



(6) 



whose uncertainty is 23 ppb. 

The last method for a that we want to men- 
tion here is based on the comparison of the mea- 
surement [llj of the fine structure of helium with 
QED predictions [12], having at present a 31 ppb 
uncertainty, 

cT^He) = 137.036 001 100(3900) th (1600) OX p. (7) 

In addition to the two methods of determin- 
ing the electron-to-proton mass ratio that were 
already mentioned, there are also ongoing efforts 
employing simple molecules. Rotational and vi- 
brational excitations have energies proportional 
to m e /M and y/m e /M, respectively. Here M 
denotes the mass of the nuclei. Thus, measure- 
ments of rovibrational spectra of systems such as 
ions H^ and HD + directly access m e /m p [TB"] . 
This approach depends of course on the accurate 



theoretical description of the ion. The current 
status of experimental and theoretical studies is 
described in |14ll5j . 

2. Bound electron g factor 

2.1. Binding effects in g — 2 

We now focus on one system: an electron 
bound in a hydrogen-like ion. We discuss its theo- 
retical description, explain how the ion is used to 
determine the electron-to-proton mass ratio, and 
discuss prospects for improvements of the theo- 
retical predictions. 
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Figure 1. Sample diagrams for the g factor of 
a bound electron at order a(Za) 5 . Solid and 
wavy lines denote electrons and photons. The 
dotted lines denote the interaction with the nu- 
cleus which is described in Sec. [31 The interaction 
vertex of the external photon is denoted with the 
symbol ®. 



Consider an ion consisting of a nucleus with 
zero spin and a single electron. We will be inter- 
ested here in the ground state of this ion. Since 
the nucleus has no magnetic moment, the mag- 
netic interaction of this system is very similar to 
that of a free electron, modulo small binding ef- 
fects. We define the bound-electron gyromagnetic 
ratio g by the correction to the energy of the ion, 
linear in the magnetic field B. In the ground 
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state, the whole magnetic moment fi of the ion is 
due to the electron's spin s, so the energy shift is 



Radial integrals of this type are 



5E = -{n-B) 



2m, 



(s-B), 



(8) 



where — e is the electron's charge. The expecta- 
tion value is to be calculated with the wave func- 
tion of the ion. 

In the lowest order, neglecting all effects ex- 
cept the Coulomb field, g was calculated [16] al- 
ready in 1928, in one of the first applications of 
the Dirac equation [T7]. To illustrate the influ- 
ence of binding effects, we reproduce here that 
result, treating the interaction with the external 
magnetic field as a perturbation. A uniform field 
is described by the vector potential A = ^B x r, 
and the correction to the energy is 



5E (Coulomb) — e cPripA ■ jip 



(9) 



where tp is the solution of the Dirac equation with 
the Coulomb potential V (r) = — — . Here Z is 
the number of protons in the ion's nucleus and 
r describes the position of the electron relative 
to the nucleus (we assume the nucleus to be very 
heavy and motionless). The ground state wave 
function can be written in the form 1181 



ip(r) = f (r) (1 + Z77 5 £ • r) v , 
f(r) = Nr- Za ' f cx-p(-Zam e r), 



(10) 



7 



Zct 



Za 



where v is a spinor such that ivv = 1 and N is a 
normalization constant such that 



(1 + 7 2 ) / d 3 rf 



(11) 



N approximately equals the value of the non- 
relativistic Coulomb wave function at the origin 

N 2 ~ ■ am< "' . In the standard (Dirac) repre- 

a ) and 7 = { 1 
With this wave function, we find 



sentation, S 



5E=—e 1 {v*B-av) \ 
3 Jo 



drr 6 f 2 (r) 



(12) 



drr n f z (r) 

{2Zam e f~ n r (to + 1 - 2Za 7 ) 



4tt (1 + 7 2 ) T (3 - 2Zaj) 



leading to 



g (Coulomb) = - I 1 + 2\j 1 - {Za) 
the result first obtained by Breit [IB] . 



(13) 



(14) 



2.2. Mixed binding-selfinteraction effects 

In view of recent precise measurements with 
hydrogen-like ions, a variety of corrections have 
recently been studied, including electron self- 
interaction loops [19I20J . vacuum polarization 
J21I22I23J . and recoil corrections [24125] . Some 
effects have been determined to all orders in Za 
[26127128129] while for others a few terms of the 
expansion around small Za have been calculated. 

Here we take the nucleus to be infinitely heavy 
and neglect small hadronic and weak effects. The 
theoretical prediction for g in this limit can be 
organized as a double series expansion in a/it 
(describing electron selfinteraction) and Za, de- 
scribing interactions between the electron and the 
nucleus, i.e. the binding effects. 

If we neglect selfinteractions, g is given by 
Breit's result ([i"4]) . On the other hand, if we ne- 
glect the binding, g is the same as for a free elec- 
tron, g = 2 + 2a, with the anomaly a described 
by eq. Q. 

The interplay of the two types of effects occurs 
first at the order O (^(Za) 2 ). Interestingly, at 
this order in Za, the correction is universal to all 
orders in a. It is found by considering the mag- 
netic interaction of a bound particle, taking into 
account its anomalous magnetic moment. The 
result is [30 31 Sg — 2a 1 1 + g ) , where a is 

given in ([T]) . 

There are no effects linear or cubic in Za, 
but beginning with (Za) 4 " all higher powers are 
present. In fact, in the approximation of a single- 
photon electron selfinteraction, it is possible to 
use the exact form of the electron propagator in 
the Coulomb field and numerically determine ef- 
fects of O (%(Za) n ) to all orders n [26127128] . 
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It is still very interesting to compute analyt- 
ically the coefficients of the Za expansion to 
check the numerical evaluation. Such a calcu- 
lation was performed for the terms —(Za) and 
f (ZafinZa only in 2004 [19], which illustrates 
the difficulty of performing even one-loop calcu- 
lations for bound particles. 

The situation becomes even harder when two 
selfinteraction loops are included. In this case 
no numerical all-order study has been performed. 
The terms (f f (Za) 4 and (f f (ZaflnZa are 
already known [20] . Together with measurements 
with carbon ions [32] the resulting theoretical pre- 
diction provides the best value of the electron-to- 
proton mass ratio, 



— ^ = 1836.152 672 9(10) 



(15) 



where we have used the CODATA recommended 
value of the proton mass in atomic units, m p = 
1.00727646677(10)u, which contributes an order 
of magnitude less to the uncertainty in ([i~5]) than 
the electron. We see that the error in (fT5|) is four 
times smaller than in the determination using a 
free electron, eq. ([3]). 

The uncertainty in (|15[) is dominated by the 
experiment; the theory is estimated to contribute 
only about three per cent of the total error. The 
reliability of the theoretical prediction is based on 
the numerical knowledge of higher-order binding 
effects in the one-loop selfinteraction. 

It is thus very important to independently 
check the numerical studies by computing fur- 
ther terms in the Za expansion. In particular, 
the coefficient of the contribution — (Za) 5 to g 
(for the principal quantum number n = 1, that 
is the ground state) is estimated to be H\(Z = 
0) = 23.15(10) [33], in the notation of [26]. 

The two- loop effect (^) (Za) 5 is estimated to 
have the coefficient between —118 and —75 j23j . 

In order to verify these predictions, we have 
undertaken an explicit determination of the coef- 
ficients of (Za) 5 terms at one and two-loop self- 
interaction level. In the next section we present 
some details of our method with examples of sim- 
pler observables, not involving an external mag- 
netic field: the Lamb shift in hydrogen-like ions, 
and the hyperfine splitting (HFS) in muonium. 



3. Lamb shift and hyperfine splitting 

In this section we describe our calculations of 
the radiative-nonrecoil corrections to the Lamb 
shift [31] and the hyperfine splitting [35] of 
hydrogen-like atoms. In both cases the shift of 
the energy levels is given by a delta function po- 
tential, which affects only S states, 



5E = -M\iJj(0)\ 2 = -M 



(Zanf 



(16) 



-0(0) is the wave function of a bound S state with 
principal quantum number n and reduced mass 
\i at the origin. The amplitude Ai is determined 
from the corresponding Feynman diagrams. Sam- 
ple diagrams are shown in Fig. [5] 
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Figure 2. Sample diagrams for the Lamb shift 
and hyperfine splitting at order a 2 (Za) 5 . Solid 
and wavy lines denote electrons and photons. The 
dotted lines denote the interaction with the nu- 
cleus (cf. Fig OH). 



The interaction of the electron with the nucleus 
is depicted in Fig. [3J In our calculation we con- 
struct an expansion in the ratio of electron and 
nucleus mass, using the method of expansion by 
regions [36]. We are only interested in the lead- 
ing term, which is given by the region where all 
loop momenta are of the order of the electron 
mass. Expanding the nucleus propagators wc find 
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Figure 3. The interaction with the nucleus (dou- 
ble line) consists of the sum of direct and crossed 
photon exchange. 



that to leading order the only difference between 
the diagrams with direct and crossed photon ex- 
change is an overall sign and the sign of the ie 
prescription of the latter. We have 



1 



1 



(q + k) 2 - M 2 + ie (q - k) 2 - M 2 + ie 
1 1 



2q ■ k + ie 2q ■ k — ie 
-2in8(2q-k) + ... , 



(17) 



where the ellipses stand for higher order terms, q 
and M are the momentum and mass of the nu- 
cleus with q 2 — M 2 , and k is the loop momen- 
tum. In the rest frame of the nucleus, the factor 
S(2q ■ k) = S(ko)/2M ensures that the photons 
exchanged between the electron and the nucleus 
carry no energy. The integration over the loop 
denoted with the dotted line in the Feynman dia- 
grams in the Figures is only over the spatial com- 
ponents k. 

Our calculation is automated to a large ex- 
tent. We use the program QGRAF [57] to gen- 
erate all Feynman diagrams. These are turned 
into FORM-readable [55] expressions by q2e and 
exp [39140] . The program MATAD3 [41] is used 
to perform the traces and express the ampli- 
tude M in terms of scalar integrals, using cus- 
tom made routines. Finally, we use integration- 
by-parts identities [42143] to express all scalar in- 
tegrals in terms of 32 so-called master integrals. 
This is achieved with the help of the program 
fire [33] , which implements the so-called Laporta 
algorithm [45146) . Results for all master inte- 
grals and details of their calculation are given in 
Ref. EH- 



4. Conclusions 

Our results for the Lamb shift and hypcrfine 
splitting at order a 2 (Za) 5 are 



SE 



Lamb 



-6.86100(4) 



a 2 (Zaf ( fj, 



irn 3 V m 



SE RFS = 0.77099(2) 



a 2 {Za) 



E F 



m e 

(18) 
(19) 



respectively. Ep = (8/z 3 (Za) 4 3Ar)/(6m e M) is the 
Fermi energy, g^ is the gyromagnetic factor of 
the nucleus. Our results improve the precision 
of the previous results in [47j by more than an 
order of magnitude. Results for individual dia- 
grams, including new analytical results, are given 
in Refs. [34135] . 

Our new results agree very well with [47j , which 
gives us confidence in the values of master inte- 
grals we have derived. The same set of 32 inte- 
grals is sufficient to compute the bound electron 
g factor to the desired accuracy, (— ) (Za) 5 . The 
only remaining obstacles are the larger number of 
diagrams and the necessity of determining correc- 
tions to the wave function of the ion. 

As the first step, we have found the gauge in- 
variant set of one-loop (vertex) corrections to the 
electron interaction with the magnetic field, an 
ingredient of the coefficient of — (Za) 5 , providing 
approximately a quarter of the estimated value of 
23.15(10) [33 . Work on the remaining one and 
two-loop diagrams is in progress. 
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